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Abstract 

The doping-induced metal-insulator transition in two-chain systems of 
correlated fermions is studied using a solvable limit of the t — J model and 
the fact that various strong- and weak-coupling limits of the two-chain 
model are in the same phase, i.e. have the same low-energy properties. 
It is shown that the Luttinger-liquid parameter K p takes the universal 
value unity as the insulating state (half-filling) is approached, implying 
dominant d-type superconducting fluctuations, independently of the in- 
teraction strength. The crossover to insulating behavior of correlations as 
the transition is approached is discussed. 
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Models of two parallel chains of correlated fermions are of interest for the 
understanding of the physical properties of a number of systems: (i) two-leg spin 
ladder systems of the type Sr x Cai4_a;Cu2404i become conducting and even su- 
perconducting under pressure, |Q even though they are insulating under ambient 
conditions; (ii) the electronic structure of certain types of carbon nanotubesJJ 
is described by very similar models; |4| (iii) two-channel quantum wires may also 
show interesting interaction effects. The two-chain case, being much easier to 
treat by controlled analytical and numerical methods than genuinely two- or 
three-dimensional models, is also helpful in understanding long-standing ques- 
tions about the existence of superconductivity in models of correlated fermions. 

Theoretical work on the two-chain model has either considered situations 
well away from half-filling, 0, [], [7| or concerned the half-filled (one electron per 
site) case where umklapp scattering leads to a Mott insulator state with a spin 
gap.|| |S], |7|, [TtJ However, very little work JTD|, |TT[ concerned the close vicinity 



of half-filling, a situation of considerable interest, in particular in view of the 
physics of doped spin ladders. In the present paper I wish to show that in that 
situation for quite general interaction strengths the Luttinger liquid parameter K p 
becomes unity (the same as for noninteracting fermions) and that consequently 
one can expect a "d-wave" superconducting state. This is to be contrasted with 
the single-chain case where close to half-filling one finds the universal strong 
coupling value K p = 1/2 and dominant antiferromagnetic correlations. fT2|, |13| , |H 



Numerical work has shown strong indications of superconductivity in Hubbard 
and t — J ladders, both calculating correlation functions directly |T5], [IB], [T7|, [TS], [05 



and by exploiting the periodicity of the ground state energy in the presence of 
a magnetic flux. EO] Nevertheless, numerical results are very hard to interpret 



in the vicinity of half-filling and numerical work has thus not helped much in 
understanding the doping-induced metal-insulator transition. 

I will now summarize the current understanding of the properties of two chains 
of interacting fermions coupled by interchain hopping and well away from half- 
filling. In the case of weak interactions a combination of perturbative renormal- 
ization group calculations and bosonization 0, |6|, [7j has shown the existence of a 
gap in the spin excitation spectrum and that superconducting correlations of d 
type dominate. More precisely, define the superconducting order parameter as 

A r = —j= (ai| r a 2 | r — O-lir^r) ■ (1) 

Here r labels sites along the chains, and ai sr is the fermion annihilation operator 
on site r, chain i, and with spin projection s. It is appropriate to call this 
order parameter "d-wave" because in Fourier space components with transverse 
wavevector and n have opposite sign. Correlations of this order parameter 
decay slowly with distance: 

(A+A ) « r"^, (2) 
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with rjscd = l/(2ifp) and K p a number close to unity. Powerlaw decay of corre- 
lations also exists for the 4k p component of the particle density: 



(n r n ) ps cos(4/c F r)r VAk F , (3) 
and one has the scaling relation 

VscdV4k F = 1 (4) 

(one has kp = ixn/2, where n is the average number of fermions per site). All 
other correlation functions, in particular those representing magnetic ordering, 
decay exponentially with distance. 

Remarkably, analogous results can also be obtained for some strong-coupling 
cases: first, if correlations within a single chain are so strong that K p < 1/3, 
renormalization-group generated interchain interactions dominate over the single 
particle hopping, and the resulting problem can be solved, leading to the same 
powerlaws (0,0) and the same scaling law (Q) as in the weak coupling case. How- 
ever, because now K p < 1/3, the Akp CDW correlations dominate over d-wave 
superconductivity. Secondly, a particular limit of the two-chain t — J model, to 
be discussed in more detail below, also leads to the same powerlaws and the same 
scaling relation. The natural conclusion from these findings in three different lim- 
its is that the correlated two-chain model is in the same phase, characterized by 
eqs. ([|JJ ; [|j in a large region of interaction strength, both for weak and for strong 
correlation. 

I now wish to determine the value of the Luttinger liquid parameter K p in 
the vicinity of the doping-induced metal-insulator transition, i.e. for n — > 1 but 
n^l. The direct use of the renormalization group approach is highly impractical 
because in this case umklapp interactions play an important role at high energies, 
but drop out at low energies. At intermediate energies, one then passes through 
a strong-coupling region which is impossible to treat systematically. The u K p < 
1/3" approach is equally impossible because for n — > 1 in a single chain one 
universally has K p — > 1/2, [|1^, [L3|, [L4| i.e. one drops outside the validity range of 
this approach. However, the t — J model gives us the possibility to reach some 
exact and general conclusions. To be specific I consider the spatially anisotropic 
u t — J — J± model" [|HJ with Hamiltonian 

H = -tJ2( a tsr a isr+i + H.c.) - t± ^2(af sr a 2 sr + H.c.) 

i,r r 

(5) 

i,r r 

where S ir = af sr a is ' r (cr) ss i/2 is the spin operator on site (i,r) and the usual no- 
double-occupancy constraint is imposed. J and J± are the exchange constants 
along and perpendicular to the ladder, and t and t± the corresponding hopping 
integrals. 
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Figure 1: Low-energy Hilbert space of the t — J — J± model in the limit J± ^> 
t,t±,J: rungs are either doubly occupied and in a singlet state (full dots), or 
both sites are empty (shaded ellipses). 



In the limit where J± ^> J, t, t± analytical progress can be made:[jn| then the 



low-energy Hilbert space consists of rungs of the ladder where either both sites 
are occupied and form a singlet or both sites are empty (figJT|). Singly occupied 
sites or rungs with a triplet lie higher by an energy of order J±. It is then 
convenient to consider the state where all rungs are occupied by singlets as the 
physical vacuum and creation of an empty rung as creation of a boson on site r, 
with associated boson creation operator 6+ = A r . In second order perturbation 
theory in t and J one then obtains an effective Hamiltonian for the low-energy 
Hilbert space: 

H eS = -t e g ^2(b+b r+1 + H.c.) + Vcs n r n r+ i , (6) 

r r 

where t eff = 8t 2 /3J±, V c s = (16t 2 /3 — 3J 2 /8)/Jj_, n r = b+b r , and now there is 
the "hardcore constraint" n r < 1. Because of that constraint, one can straight- 
forwardly transform the model into a spin-1/2 language by setting = b£ and 



S* = n r — 1/2, to recover the well-known spin-1/2 XXZ spin chain. [[H], |22| In 
Haldane's paper the interesting Luttinger liquid parameter K p is determined. 
In particular, if the original fermionic model is close to the insulating state at 
half-filling, n — > 1, corresponding to a very dilute hardcore boson gas, one has 
K p — 1 + 0(1 — n), independent of the value of the interaction, i.e. K p takes a 
universal value when the metal-insulator transition is approached. Moreover, I 
have argued above that both the weak-interaction limit of the two-chain model 
and the strongly interacting chains weakly coupled by interchain hopping are in 
the same phase as the t — J — J± model. One thus concludes that upon approach- 
ing the doping induced metal-insulator transition, n — ► 1, the Luttinger liquid 
parameter K p takes the universal value unity, independently on whether correla- 
tions are weak or strong. Via the relations ([§|3],H) this then implies that d-type 
superconducting correlations are strongly dominant close to the metal-insulator 
transition. 

It is interesting to understand how the powerlaws (0,0) connect to the behav- 
ior expected for the insulating case. Consider first the Akp CDW correlations: at 
distances shorter than the average spacing between bosons one expects the be- 
havior typical of an insulator, i.e. a constant. A reasonable form for the crossover 



between short-distance insulating behavior and the asymptotics of eq.@ then 
isfl23f| 



(n r n ) « cos(4k F r)(r + r 2 )~ Kp , (7) 

where £ oc 1/(1 — n) is a length proportional to the average distance between 
bosons. On the other hand, the asymptotics @ of the pairing correlations are 
expected to be valid only at distances larger than £, with exponential decay at 
shorter distances. This in particular implies that the amplitude of the divergence 
of the pairing susceptibility (xscd oc max(cj, q, T)~ 3 / 2 as K p — > 1) vanishes as 
n — > 1. 

Recently, Konik et al. ||ll|| have investigated the vicinity of the metal-insulator 
transition in the weak-coupling limit, starting from the 5*0(8) model of the in- 
sulating state [TIJ and using the exact integrability of that model. They arrive 
at the same conclusion, K p — 1 + 0(1 — n). This provides a confirmation of the 
continuity between weak and strong correlation conjectured before. 

I finally comment on some properties of the model (|5[) for general fermion 
density. First, in the very dilute limit, corresponding to a dense bosonic model 
(P), from Haldane's results one again has K p = 1 + 0(n). More interestingly, 
for a quarter-filled fermionic model (one boson per two sites), there is a metal- 
insulator transition into an insulating CDW state with period 4k p when V^r > 
2i e ff. In the XXZ spin chain, this transition is due to umklapp processes which 
in the original fermionic language of eq. (|5]) correspond to four-fermion umklapp 



processes, as first discussed in ref.[23], and which can become relevant only for 
strongly repulsive interactions. In the conducting state (V c q < 2t eff ) one has 
K p > 1/2, i.e. superconductivity dominates. On the other hand, for V e s > 2t eff , 
when n —> l/2 ± one has K p — > 1/4, i.e. 4kp CDW fluctuations dominate. The 
t— J — J± model thus provides an interesting example for the transition between 
dominant superconducting and CDW fluctuations. Finally, we remark that in the 
case of strongly attractive interactions, V c s < —2t e s, one has phase separation 
for any filling of the band. 

In conclusion, based on a solvable limit of the t — J ladder model and a con- 
tinuity conjecture between strong and weak correlation, I have shown here that 
close to the doping-induced metal-insulator transition the Luttinger liquid pa- 
rameter K p of the two-chain model takes the universal value unity, corresponding 
to dominant d-type superconducting correlations. This contrasts strongly with 
the single-chain case, where K p — > 1/2 and antiferromagnetism dominates. It 
should however be pointed out that superconductivity is rather easily destroyed 
by disorder, |23[] and that therefore in any real system the existence or not of 



superconductivity will depend crucially on the interplay between disorder and 
interladder coupling which stabilizes superconductivity. 



5 



References 

M. Uehara et al, J. Phys. Soc. Jpn. 65, 2764 (1996). 
[2] E. Dagotto and T. M. Rice, Science 271, 618 (1996). 
[3] T. Ebbesen, Phys. Today 49, 26 (1996). 
[4] C. Kane and E. J. Mele, |cond-mat/9608146| . 
[5] M. Fabrizio, Phys. Rev. B 48, 15838 (1993). 
[6] H. J. Schulz, Phys. Rev. B 53, R2959 (1996). 
[7] L. Balents and M. P. A. Fisher, Phys. Rev. B 53, 12133 (1996). 

H. J. Schulz, Phys. Rev. B 34, 6372 (1986). 
[9] S. P. Strong and A. J. Millis, Phys. Rev. Lett. 69, 2419 (1992). 



[10] H.-H. Lin, L. Balents, and M. P. A. Fisher, preprint |cond-mat/9801285 
(unpublished) . 



[11] R. Konik, F. Lesage, A. W. W. Ludwig, and H. Saleur, preprint |cond- 
| mat/ 9806334| (unpublished). 

[12] H. J. Schulz, Phys. Rev. Lett. 64, 2831 (1990). 

[13] H. J. Schulz, Phys. Rev. B 22, 5274 (1980). 

[14] T. Giamarchi, Phys. Rev. B 44, 2905 (1991). 

[15] E. Dagotto, J. Riera, and D. J. Scalapino, Phys. Rev. B 45, 5744 (1992). 

[16] R. M. Noack, S. R. White, and D. J. Scalapino, Phys. Rev. Lett. 73, 882 
(1994). 

[17] R. M. Noack, S. R. White, and D. J. Scalapino, Physica C 270, 281 (1996), 
bond-mat /960 1047 . 

[18] M. Troyer, H. Tsunetsugu, and T. M. Rice, Phys. Rev. B 53, 251 (1996). 



[19] K. Kuroki, T. Kimura, and H. Aoki, |cond-mat/9610038| . 
[20] C. A. Hayward et al, Phys. Rev. Lett. 75, 926 (1995). 
[21] C. N. Yang and C. P. Yang, Phys. Rev. 150, 327 (1966). 
[22] F. D. M. Haldane, Phys. Rev. Lett. 45, 1358 (1980). 



6 



[23] H. J. Schulz, in Strongly Correlated Electronic Materials: The Los Alamos 
Symposium 1993, edited by K. S. Bedell et al. (Addison-Wesley, Reading, 
Massachusetts, 1994), p. 187, |cond-mat/9412036| . 

[24] E. Orignac and T. Giamarchi, Phys. Rev. B 53, 10453 (1996). 



7 



